The Effect of Neutral Atoms on Capillary Discharge Z-Pinch 
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We study the effect of neutral atoms on the dynamics of a capillary discharge Z-pinch, in a regime 
for which a large soft-x-ray amplification has been demonstrated. We extended the commonly used 
one-fluid magneto-hydrodynamics (MHD) model by separating out the neutral atoms as a second 
fluid. Numerical calculations using this extended model yield new predictions for the dynamics of 
the pinch collapse, and better agreement with known measured data. 
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Z-pinch collapse has been extensively studied since the 
late 50s, being a simple and effective way of producing 
hot and dense plasma. In this process, an electric cur- 
rent flowing through a plasma column, interacts with its 
self magnetic field, and the resulting force contracts the 
plasma column in the radial direction. Today Z-pinch 
plasma is widely used for various applications such as 
high power radiation sources and neutron sources 
An exciting new application of Z-pinch plasma was re- 
cently demonstrated by Rocca et. al. |3ra]. In this work, 
large amplification of soft-x-ray light was obtained in Ne- 
like Ar and S plasma, created by a fast (~ 40 ns) Z-pinch 
discharge inside a capillary. Compared with the alterna- 
tive approach of laser driven excitation the capil- 
lary discharge has the advantage of allowing for compact 
(table-top), efficient and simpler soft-x-ray lasers. 

In this paper we study the role of neutral atoms in the 
dynamics of a capillary discharge Z-pinch, in the regime 
for which soft-x-ray amplification was demonstrated. 
The commonly used one-fluid magneto-hydrodynamics 
(MHD) model assumes that all the particles in the plasma 
are charged, and drift together. We, however, show that 
for the case discussed here, large portions of the plasma 
contain an appreciable amount of neutral atoms. Since 
these are not affected by the electro-magnetic forces, but 
only by the much weaker mechanical forces, they flow 
with much smaller velocities than the ions and the elec- 
trons. To account for this effect, we extend the one- 
fluid MHD model by introducing a separate fluid for the 
neutral atoms (in addition to the standard electrons-ions 
fluid). Results of calculations using this extended model 
give new predictions for the dynamics of the pinch col- 
lapse, with some features in better resemblance with the 
measured data. This confirms our previously reported 
estimates [||. 

We start with the standard one-fluid two-temperature 
MHD model, commonly used for numerical calculations 
of Z-pinch processes p-[l3|. It considers hydrodynamic 
flow including shock waves, heat conduction, heat ex- 
change (between ions and electrons), magnetic field dy- 
namics, magnetic forces, Ohmic heating, radiative cool- 
ing and ionization. We use a simple ionization model, 
and assume a quasi steady state, taking into account col- 



lisional ionization, and 2-Body and 3-Body recombina- 
tion. Since the plasma is assumed to be optically thin, 
ionization and excitation by radiation are neglected. The 
latter assumption should hold at least to the end of the 
collapse. This model is incorporated into our numerical 
code, SIMBA, where the equations of motion of the sys- 
tem (see [PH14[ ) are solved in a Lagrangean mesh pB), 
assuming one-dimensional axial symmetry. 

Shown to be remarkably stable fl6f| , and having a high 
length-to-diameter ratio (of 50-500), the capillary dis- 
charge Z-pinch experiment is naturally described in the 
framework of this 1-D MHD model. Previously reported 
works [jl0| , have indicated that taking into account abla- 
tion of the plastic capillary wall is necessary for the cor- 
rect description of the pinch dynamics. According to this 
the calculation should thus be extended to include a nar- 
row region of the plastic capillary wall. However, it was 
also shown in |l0| that even with this effect taken into 
account, good agreement with the measured data still 
requires some major artificial adjustments of the plasma 
transport coefficients. We have repeated these calcula- 
tions using the same one-fluid MHD model, and found 
them to agree with the reported results. In particular, 
we also find that the measured data is reproduced by one- 
fluid MHD calculations only when artificial adjustments 
are introduced, as demonstrated in Fig. (|l|). The figure 
displays the calculated radius of the collapsing Ar plasma 
as a function of time in a capillary discharge Z-pinch. The 
parameters of the calculations are those used for soft-x- 
ray amplification experiments initial Ar density 
of po = 1.7- 10 _6 g./cm 3 or Uq ~ 2.5T0 16 atoms/cm 3 , ini- 
tial temperature of To « 0.5eV, and a maximum current 
of 39kA, with its peak at t=32ns p7| . The figure also 
presents some measured data, of the radius of soft-x-ray 
source, as a function of time, taken from [Q. Since the 
radii of the soft-x-ray source and that of the collapsing Ar 
plasma are related, it is clear that there are disagreements 
between the calculated and measured data: For example, 
the calculated pinch peak is about 10ns earlier than the 
measured one. It is shown in Fig. M) that multiplying 
the classical electrical conductivity p"q] by a factor of 15, 
results in a good agreement with the measured instant of 
the pinch peak, however at the same time it also spoils 
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FIG. 1. Z-Pinch of Ar plasma inside a plastic capillary. 
Thin line: Calculated outer boundary of Ar plasma, assuming 
classical electrical conductivity. Thick line: Calculated outer 
boundary of Ar plasma, with a factor of 15 on the classical 
electrical conductivity. Dots: Measured radius of soft-x-ray 
source [hoi. 



the agreement with measured collapse velocity. We no- 
tice that both calculations do not properly reproduce the 
initial stages of the collapse, which is delayed by about 
10-15ns. According to reproducing the whole stages 
of the measured collapse requires more artificial adjust- 
ments in the plasma transport parameters, up to 20-40 
times their classical values. This need for artificial adjust- 
ments of plasma parameters in one-dimensional one-fluid 
MHD calculations cannot be explained by two- or three- 
dimensional effects in the modeled experiment: The work 
of Bender et. al. jl6| has proven a perfect azimuthal 
(</>-direction) symmetry in this same capillary discharge 
Z-pinch, and the demonstrated amplification gain 
indicates a very good Z-direction symmetry. 

In order to better understand the dynamics of the 
pinch collapse, we have focused our study on the impor- 
tance and the role of neutral atoms in this process. The 
one-fluid MHD model assumes that the plasma consists 
of two components: electrons and effective single- type 
ions, with their charge being the average charge of all 
the differently charged ions in the plasma, including the 
neutral atoms. In addition, these two components are 
assumed to flow together, as a single fluid. This assump- 
tions are reasonable for regimes for which at least one of 
the two following conditions is fulfilled: (i) All the atoms 
in the plasma are ionized, or, (it) The neutral atoms are 
strongly coupled to the charged particles, and hence fol- 
low them in the same single fluid. 

Fig. (||) presents the percentage of neutral atoms 
as a function of electron temperature in Argon plasma, 
based on our ionization model. According to this fig- 
ure, a plasma of electron temperature lower than 2-3 eV 
contains an appreciable amount of neutrals. In Carbon 



plasma, which is a typical representative of the ablated 
capillary wall, the picture is similar. Our MHD calcula- 
tions show that the Ar plasma starts to heat up above 
2-3 eV only 5 ns after the beginning of the pinch, and its 
central region stays below this temperature for the next 
25 ns Major portions of the plastic wall plasma re- 
main below 2-3eV even after the pinch collapses at the 
axis. The percentage of neutral atoms in the plasma is 
hence far from being negligible. We thus conclude that 
condition (i) does not hold. The plasma contains three 
different components: electrons, ions, and neutral atoms. 
We now turn to check whether or not condition (ii) is sat- 
isfied, by examining the couplings between these different 
ingredients of the plasma. The electrons and ions, being 
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FIG. 2. Percentage of neutral atoms in Ar plasma, as a 
function of electron temperature. Each line corresponds to a 
different plasma density. 

charged particles, are coupled through Coulomb forces. 
A measure of the strength of this coupling is given by 
the plasma frequency, up. For the case discussed here, 
l/ujp « 10~ 5 — 10~ 3 ns, which is negligible compared to 
the typical pinch collapse times of Tp^^ ~ 40ns. This 
means that the coupling between the electrons and ions 
is very strong, and that they practically drift together, as 
a single fluid. The neutral atoms, however, are coupled 
to the charged particles only by collisions, and may thus 
flow separately, as a second fluid. We therefore assume 
two fluids, one of charged- species (electrons and ions) and 
the other of neutral- species (atoms), with flow velocities 
Ui and u a respectively. The collisional momentum trans- 
fer between these two fluids is evaluated assuming a hard 
spheres approximation: We regard the two fluids as two 
clouds of hard spheres, drifting through one another. In 
that case, the collision frequency per unit volume equals: 

C" = a-r 2 a -n a -ni-\u a - Ui \ (1) 

and the collisional momentum transfer rate, per unit vol- 
ume, is thus 



t — ct • r 



\u a - Ui\ (u a - m) , (2) 



where a is a coefficient of about 27r. Here r stands for 
the particle radius, m for its mass, and n stands for the 
number density. The indices a,i denote atoms and ions 
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respectively. Later on we will use the index e for elec- 
trons. 

The force in Eq. (||), F°f u , depends quadratically on 
the velocity difference between the charged-species and 
the neutral-species fluids. This coupling thus restrains 
the separation between the two fluids. Taking reasonable 
densities of n a sw 10 16 ,rti « fO 15 (10% ionization), and an 
appreciable velocity difference of \u a — Ui\ « 10 6 cm/s, we 
get for Ar plasma a collisional coupling term of the order 
of 10 6 dyn/cm 3 . This is 2-3 orders-of-magnitude less then 
the estimated magnetic ("j* x ~B*/c) and hydrodynamic 
(VP) forces. 

We conclude that in the regime discussed here, both 
of the above conditions for the validity of the one-fluid 
MHD fail to be satisfied. The two fluids are indeed ex- 
pected to flow separately. However, they exchange mass, 
momentum and energy due to exchange of particles (by 
ionization and recombination) and due to atoms-ions col- 
lisions. By S a (r,t) we denote the mass sink (per unit 
volume, per unit time) in the neutral-species fluid due to 
ionization of neutral atoms (S a > 0). S a plays a role of 
a source in the charged-species fluid. Similarly, Si(r,t) 
denotes the mass sink in the charged-species fluid, due to 
recombination of ions +1 (Si > 0). The total mass trans- 
fer from the neutral-species fluid into the charged-species 
fluid due to ionization and recombination is thus S a — SV 

To account for the exchange of mass, momentum and 
energy between the two fluids the standard one-fluid 
MHD for the charged-species fluid (see (li) for exam- 
ple) are amended, and new, separate equations for the 
neutral-species fluid are added. The revised mass equa- 
tion for the charged-species fluid is then (we use cylin- 
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I (ft + ft) + iBi+fAl. {rUi) = {Sa - Si ) , (3) 
at r or 

where p stands for mass density, and = + u ■ V is 
the comoving derivative. The separate mass equation for 
the neutral-species fluid is then: 

J t (P«) + ^^(ru i ) = -(S a -S i ) (4) 

The revised momentum equation for the charged species 
fluid is: 



, \ d d . , 

(ft + ft) -77^ = -5- [Pe + Pi, 
at or 
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+F a c f + S a (u a - m) 



(5) 



where P stands for pressure, j for current density and 
B for magnetic field. F^° u is the collisional momen- 
tum exchange between the neutral-species fluid and the 
charged-species fluid, given in Eq. (|^). The momentum 
equation of the neutral-species fluid should then be: 



Pa dt Ua 



Similarly, the one-fluid MHD ion-energy equation [p[-^3[) 
is also properly amended, and separate atom-energy 
equation for the neutral-species fluid is introduced. Colli- 
sions between the two fluids, as well as particles exchange 
due to ionization and recombination are considered in 
these equations in the same manner as in the mass and 
momentum equations. The MHD electron-energy equa- 
tion is left unchanged. 
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FIG. 3. Calculated outer boundary of Ar plasma during 
a Z-pinch capillary discharge. Thin line: Neutral atoms ne- 
glected (standard one- fluid MHD). Thick line: Neutral atoms 
included (extended model) . 

These equations were incorporated into our SIMBA 
code. For simplicity, and in order to emphasize the ef- 
fect introduced by separating the neutral atoms from the 
charged-species fluid, we assume, in the following calcula- 
tions, that the capillary wall is also made of Argon. The 
other pinch parameters are left unchanged, however we 
now use the classical transport coefficients |Q, without 
any artificial adjustments. Fig. (||) shows the effect of 
the neutral-species fluid on the calculated outer bound- 
ary of the collapsing Ar plasma. It is clearly indicated 
that the effect of the neutral component in the capillary 
discharge Z-pinch is not negligible. When the neutral- 
species fluid is included the collapse seems to be delayed, 
however after it starts it is more rapid. This trend seems 
to better resemble the data presented in Fig. ([[]), where 
it was shown that compared to one-fluid MHD calcula- 
tions the measured collapse is delayed, and after it starts 
the collapse rate is much higher. 

We have also examined the effect of neutral atoms on 
the electron density distribution during the pinch. In 
Fig. ([|), the calculated spatial distribution of electron 
density at time=25ns is plotted, with and without the 
neutral-species fluid. Both models predict a collapsing 
plasma sheath, and show some ablated material from the 
capillary wall. However, when the neutral-species fluid 
is taken into account, the collapsing plasma sheath is 
wider and less dense, compared to the predictions of the 
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FIG. 4. Calculated electron density profiles at time=25ns. 
of an Ar Z-pinch capillary discharge. Thin line: Neutral 
atoms neglected (standard one-fluid MHD). Thick line: Neu- 
tral atoms included (extended model). 

standard one-fluid MHD model. 

We like to offer a qualitative explanation for the re- 
sults presented in Fig. In the one-fluid MHD 
model, the atoms and ions are assumed to flow together 
with the electrons. The magnetic forces, which are domi- 
nant in this case, thus accelerate the whole plasma body. 
In reality, however, only the ions flow together with the 
electrons, while the neutral atoms flow separately. Since 
the plasma is initially mostly neutral, the magnetic forces 
act only on a small fraction of the total mass, which is 
then rapidly accelerated inwards. Most of the Ar stays 
outside, almost at rest. While the process evolves, more 
atoms get ionized, and join the charged-species fluid. 
This effect is seen in Fig. (||) as a delay in the collapse. 
At any given spatial and temporal point, the magnetic 
forces act on a "freshly" ionized matter, almost at rest. 
The resulting acceleration is thus more gradual, leading 
to a wider and less dense plasma sheath, as seen from 
Fig. (|). 

In conclusion, we have shown that the effect of neu- 
tral atoms on the dynamics of the capillary discharge Z- 
pinch is not negligible. We have demonstrated that sep- 
arating out the neutral atoms as a second fluid produces 
a different pinch collapse dynamics, with some features 
similar to the measured data. It is expected that the 
improved modeling of the pinch collapse dynamics will 
yield a better understanding of capillary discharge X-ray 
lasers, since the amplification gain, as well as the propa- 
gation and refraction of radiation in the lasing media are 
both dominated by the details of the plasma state. 
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